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Abstract 



We address the classical and quantum marginal problems, namely the question of 
simultaneous realizability through a common probability density in phase space of 
a given set of compatible probability distributions. We consider only distributions 
authorized by quantum mechanics, i.e. those corresponding to complete commuting 
sets of observables. For four-dimensional phase space with position variables q and 
momentum variables p, we establish the two following points: i) given four compatible 
probabilities for ((71,(72); {qi,Pz), (pi,Q2) and (pi,P2), there does not always exist a 
positive phase space density p(q,p) reproducing them as marginals; this settles a long 
standing conjecture; it is achieved by first deriving Bell-like inequalities in phase space 
which have their own theoretical and experimental interest, ii) given instead at most 
three compatible probabilities, there always exist an associated phase space density 
p(q,p); the solution is not unique and its general form is worked out. These two 
points constitute our "three marginal theorem" . 



PACS : 03.65.Ta, 03.67.-a 



I. Introduction 



In classical mechanics position and momentum can be simultaneously specified. Hence 
phase space density has a well defined meaning in classical statistical mechanics. In 
quantum theory the probability density for observing eigenvalues of a complete com- 
muting set (CCS) of observables is specific to the experimental context for measuring 
that CCS. Joint probabilities for different CCS which contain mutually noncommut- 
ing operators are not defined. For example for a 2-dimensional configuration space, 
with q,p denoting position and momentum, probability densities of anyone of the four 
CCS (91,^2)) {qi>P2), (pi,Q2) or (pi,P2) are defined, but not their joint probabilities. 
The question one may raise is : can one define such joint probabilities, e.g. a phase 
space probability density p(q,p) such that all its marginals^] coincide with the quan- 
tum mechanical probabilities for the different individual CCS? This general question 
was first raised by Martin and Roy 0. 

The Martin- Roy contextuality theorem demonstrates the impossibility of realizing 
quantum probability densities of all possible choices of the CCS of observables as 
marginals of one positive definite phase space density. For example, consider a two 
dimensional configuration space. Let coordinates q\ a , q2 a be obtained from q±,q2 by 
a rotation of arbitrary angle a, and momenta pi a , p2 a be related similarly to pi,P2 

qia\_y(qi\ fpia\ = yfpi\ (u) 



(12a J V<?2/ \P2a/ \P2 

where 

y _ f COS a Sm Q \ /j 2) 

\ — sin a cos a J 

Does there exist for every quantum state (with density operator p) a positive defi- 
nite phase space density p(q,p) such that its marginals agree with the corresponding 
quantum probabilities, i.e., 

J dp la dq 2a p(q,p) = (qia,P2a\p\qia,P2a) (1-3) 
for all a ranging from to 2tt? They answered this question in the negative by finding 



a state p for which eqs. (1.3) for all a are inconsistent with positivity of p. Since 
different a correspond to different experimental contexts, the Martin-Roy theorem 
is a new Gleason-Kochen-Specker type contextuality theorem [0]. The positivity of 
the phase space density p(q,p) is absolutely crucial for this theorem; otherwise the 
Wigner distribution function M would be a solution of (jL3|). 



Equations (L3) constitute conditions on an infinite set of marginals of p(q,p) 
(corresponding to the continuously infinite choices for a) to agree with corresponding 
quantum probability densities. Their inconsistency still leaves open the question of 
consistency of a finite number of such marginal conditions. 

In agreement with common terminology, by "marginal" of a distribution over several variables, 
we denote integrals of the distribution over a subset of its variables. 
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Indeed, the consistency of two marginal conditions where the marginals involve 
only nonintersecting sets of variables has been known for some time. Cohen and 
Zaparovanny Q constructed the most general positive p(q,p) obeying 

J dp p(q,p) = (q\p\q), J dq p(q,p) = (p\p\p). 
Their solutions generalize the obvious simple uncorrelated solution for pure states ip, 

p(q,P) = mq)\ 2 \m\ 2 , 

where tilda denotes Fourier transform. Based on generalized phase space densities 
exhibiting position momentum correlations, Roy and Singh || constructed a causal 
quantum mechanics reproducing quantum position and momentum probability den- 
sities, thus improving on De Broglie-Bohm mechanics || which only reproduced the 
quantum position probability densities. Later, going much further than the noninter- 
secting marginals of Cohen et al. 0, Roy and Singh |7| constructed a causal quantum 
mechanics based on a positive p(q,p) whose marginals reproduce the quantum prob- 
ability densities of a chain of N + 1 different CCS, e.g. 

(Qi,Q 2 ,---,Qn), (Pi,Q 2 ,---,Qn), (Pi,P 2 ,Q 3 ,---,Q n ),---(Pi,P2,---,Pn)- 

Here N is the dimension of the configuration space, and each CCS in the chain is 
obtained from the preceding one by replacing one of the position operators Qi by the 
conjugate momentum operator Pj. 

Roy and Singh proposed the following definition: a Maximally Realistic Causal 
Quantum Mechanics is a causal mechanics which simultaneously reproduces the 
quantum probability densities of the maximum number of different (mutually non- 
commuting) CCS of observables as marginals of the same positive definite phase space 
density. They also conjectured that for N dimensional configuration space this max- 
imum number is N + 1. 

A proof of this long standing conjecture is important for quantum mechanics where 
it quantifies the extent of simultaneous realizability of non commuting CCS. 

In this paper, we restrict ourselves to the case of iV = 2 degrees of freedom. The 
general case (N > 2) will be dealt with in a forthcoming paper. In Section II below, 
we first state the classical and quantum marginal problems and second, show 
that, given four classical compatible two-variable probability distributions, there does 
not always exist a positive phase space distribution reproducing them as marginals. 
In Section III, we develop a new tool, "the phase space Bell inequalities" , which are 
the phase space analogues of the standard Bell inequalities |J for a system of two spin- 
half particles. We use them in Section IV to prove the conjecture for four-dimensional 
phase space (N = 2), namely the impossibility of simultaneous realization of quan- 
tum probabilities of more than three CCS as marginals. In Section V, we explicitly 
construct the most general phase space distribution which reproduces probabilities of 
three CSS as marginals. These results, the three marginal theorem, are relevant 
for the construction of maximally realistic quantum mechanics. 

As our results are essentially new theorems for multidimensional Fourier trans- 
forms, they are also expected to be useful for classical signal and image processing |J. 



2 



The theorems of the present paper and their generalizations to arbitrary N [10] con- 
siderably advance previous results in the field, which have only dealt with noninter- 
secting sets of marginals (e.g. time and frequency). A summary of the results of this 
paper without detailed proofs is being reported separately [11|. 



II. Four marginal problem 

Let us consider a physical system with 2-dimensional configuration space. Let (qi,pi) 
and (<72jf?2) be a set of canonical variables in the corresponding phase space. We look 
for a (normalized) probability distribution p(qi,q2,Pi,p 2 ) such that 

p{qi,Q2,Pi,P2)>0, (HI) 



/ dp 1 dp 2 p(qi,q 2 ,p 1 ,p 2 ) = R(qi,q 2 ), (H.2) 

J dp 1 dq 2 p(q 1 ,q 2 ,pi,p 2 ) = S(q 1 ,p 2 ) , (II.3) 

J dq 1 dp 2 p{q 1 ,q 2 ,p 1 ,p 2 ) = T(p 1 ,q 2 ), (II.4) 

J dqidq 2 p(qi,q 2 ,pi,p 2 ) = U(pi,p 2 ) , (II.5) 

where the four marginals R(qx,q 2 ), S(qi,p 2 ), T(pi,q 2 ), and U(pi,p 2 ), the respective 
joint probabilities, are given. For consistency we must have 

R,S,T,U>0, (II.6) 

and 



dq 2 R(q 1 ,q 2 ) = J dp 2 S(q 1 ,p 2 ), 

dq 1 R(q 1 ,q 2 ) = J dpiT(px,q 2 ) . 

dqiS(qi,p 2 ) = f dpiU(pi,p 2 ) 

dq 2 T(pi,q 2 ) = J dp 2 U(pi,p 2 ) 



(II.7) 



We shall refer to the problem Given four distributions R, S, T and U, satisfying 
the consistency conditions, does there always exist a positive p(qi,q 2 ,Pi,p 2 ) with 
these distributions as marginals? as the Classical four marginal problem. 

When the system is quantum mechanical and is described by a state vector |^), 
each of the four marginals involves a pair of compatible observables and we have 

R(qi,q2) = I <<?i, ^l^lQ'i, "?2>| 2 , 

S{qi,p 2 ) = \{qi,p 2 \p\qi,P2)\ 2 , 

T(pi,q 2 ) = \(pi,q2\p\pi,q2}\ 2 , 

U(pi,p 2 ) = \(pi,p 2 \p\pi,p 2 )\ 2 • 
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In this case, the above consistency conditions are automatically satisfied. We then 
refer to the problem as the Quantum four marginal problem. A positive answer to it 
for all states p would mean that a realistic interpretation of the quantum results is 
possible (to the extent that only measurements connected to the four marginals are 
involved) . 

We shall see that the answer to both problems is negative. 

Let us first show that the classical four marginal problem does not always admit 
a solution. To this end, consider the following set of marginals 



R(qi 


Q2) 


1 

2 


[%1 


- ai)5(q 2 


- 0-2) 


+ %i 


- a' l )S(q 2 


~ a' 2 )] , 


(II.9) 


S(qi 


P2) 


1 

2 


[%1 


- ai)5(p 2 


-62) 


+ %i 


- a' 1 )S(p 2 


" b' 2 )] , 


(11.10) 


T(pi 


Q2) 


1 

2 


[S(pi 


- h)8(q 2 


- a 2 ) 


+ 5{ P i 


- b' 1 )S(q 2 


-4)] , 


(11.11) 


U(pi 


P2) 


1 

2 


[S(pi 


- h)S(p2 


~b' 2 ) 


+ S( P i 


- b' 1 )S(p 2 


- fe)] , 


(11.12) 



which obviously satisfy the consistency conditions and ( |IL7 ). They possess two 



essential features. First, their non factorized form. Second, in view of the expressions 
of R, S and T, the positions of the factors 6(p 2 — b 2 ) and 8{p 2 — b' 2 ) in the expression 
of U are not the "natural ones" . 



Eq.( p.9| ) means that the support of the distribution R in the plane (gi, q 2 ) consists 
in the two points (01,02) and (a'^a^)- As a consequence, any positive p satisfying 
flH.2p should have a support the projection of which on the plane (qi,q 2 ) would also 
consist in those two points. That is 

P = %i - a 1 )8(q 2 - a 2 )a(pi,p 2 ) + 8(qi - a[)5(q 2 - a' 2 )a (p 1 ,p 2 ) , (11.13) 

where a and a' are some positive distributions. Similarly, from eqs,( |II.3j ) to ( p.5| ) 

P = 8(q 1 -ai)5(p 2 -b 2 )P(pi,q 2 ) + 5{qi-a 1 )5(p2-b 2 )(3'(pi,q 2 ), (11.14) 
= S(pi-bi)5(q 2 -a 2 )'y(qi,p 2 )+S(pi-b' 1 )6(q 2 -a , 2 )j'(q 1 ,p 2 ), (11.15) 
= 8{pi - h)S(p 2 - b' 2 )n(qi,q 2 ) + 8(j>\ - b' 1 )5(p 2 - b 2 )rf{qi,q 2 ) . (11.16) 

According to eqs. flH.13D to (|II.15| ) 



v5(qi - a 1 )5(q 2 - a 2 )5(pi - h)5(p 2 - b 2 ) 
+v'S( qi - a' l )5(q 2 - a' 2 )5{p 1 - b' 1 )5(p 2 - b' 2 ) 



(11.17) 



with v > 0, v' > 0, (v+v' = 1). Clearly, eqs.( p,16 ) and ( [EI. 17 ) are incompatible, which 



establishes the non existence of p, and settles the classical four marginal problem. 

This however does not settle the quantum problem. Actually, the above example 
obviously cannot be strictly realized through a wave function in accordance with 



eqs.(IL8). More than that, this example is so "twisted" that, even after smoothing 



out the 5 measures in eqs.( [II.9|) to (11.12), approaching it close enough through a 



wave function appears as very difficult (if not impossible). Instead, we develop a new 
mathematical tool. 
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III. Phase space Bell inequalities 



Consider any choice of functions r(qi,q 2 ), s(qi,p 2 ), t(pi,q 2 ), and u(pi,p 2 ), obeying 
A < r(q 1 ,q 2 ) + s(q 1 ,p 2 ) +t(pi,q 2 ) + u(p 1 ,p 2 ) < B (Vqi,q 2 ,Pi,p 2 ), ( IIL1 ) 

Multiply by p(q,p), integrate over phase space and use positivity and normalization 
of p(q,p). We deduce that the (classical as well as quantum) four marginal problem 
cannot have a solution unless 

A < J dq 1 dq 2 r(q 1 ,q 2 )R(q 1 ,q 2 ) + J dqidp 2 s(q 1 ,p 2 )S(qi,p 2 ) 

+ J dp 1 dq 2 t(p 1 ,q 2 )T(pi,q 2 ) + J dp 1 dp 2 u(p 1 ,p 2 )U(p 1 ,p 2 ) < B. 

(111.2) 

Here R, S, T and U are defined by eqs. flH.cf ) in the quantum case. It turns out that 
a particularly interesting choice is 

r(qi,q 2 ) = sgniq(gi) sgaF 2 (q 2 ) , 

s(qi,P2) = sgnFi(gi) sgnG 2 (p 2 ), (HIS) 

t(pi,Q2) = sgnGi(pi) sgnF 2 (<2 2 ) , 

u(pi,p 2 ) = -sgnGi(pi) sgnG 2 (p 2 ) , 

with A = —2, B = +2 and with F\, F 2 , G\ and G 2 arbitrary non vanishing functions^. 
Then the inequalities ( |III.2j ) become a phase space analogue of the Bell inequalities 
for spin variables. 

The necessary conditions ( [II.2[ ) provide us with an alternative proof that the 



classical problem does not always admit a solution. Indeed, it is readily seen that 
they are violated for the marginals flII.9D to ( [11.12 ) and functions F's and G's such 
that 

F 1 (a 1 ),F 2 (a 2 ),G 1 (b 1 ),G 2 (b 2 )>0, 
F l (ai),F 2 (a' 2 ),G 1 (b , l ),G 2 (b' 2 )<0. 



We shall see in the next section that the necessary conditions ( III.2j ) can be violated 



also in the quantum case. There, the analogy between our correlation inequalities 
( HI-2j ) (with the choice ( [III.3D ) and Bell inequalities will become more apparent, 



especially as regards to their implications. 

IV. Solving the four marginal quantum problem 

This section is divided into four parts. In the first one, we prove the existence of 



wave functions which violate the correlation inequalities ( III. 2 ). Strictly speaking 



this already settles the problem. However, the explicit construction of such wave 

2 Note that, with this choice, the sum r + s + t + u assumes only its two extremal values A — —2 
and B = +2, which makes it in a sense optimal. 
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functions, which we present in subsections B and C, is worthwhile in that it exhibits 
the physical implications of our inequalities. In subsection D, we elaborate on the 
formal analogy with Bell inequalities. 



IV-A. Non constructive proof 

One first notices that Xi(^i) = \ [1 + sgniq(gi)] is the characteristic function of some 
set Si C K, and similarly for F 2: G\ and G2, so that eqs.( IIL^ ) read 





= (2X1 


-1)(2X2 


- 1 


s(qi,P2) 


= (2X1 


-1)(2X' 2 


- 1 


t(pi,q2) 


= (2xi 


- 1)(2X2 


- 1 



(IV.l) 

«(P1,P2) = -(2Xi-l) (2X2-l)i 
where Xi stands for Xi(%) an d Xi f° r Xi(Pi)i = 1 5 2). Inequalities ( III.1| ) then become 

< V < 1 , (IV.2) 
and in fact V(l — V)=0 (see footnote |2|), where V(qi, q2,Pi,P2) is given by 

V = Xi + X2 + X1X2 - X1X2 - X1X2 - X1X2 • (IV.3) 
Let us define a corresponding quantum operator V by 

V = Xi + X2 + X1X2 - X1X2 - X1X2 - X1X2 , (IV.4) 

where 



XI = / d qi \qi)(qi\ ® I2 , 

JSi 

X2 = ll <8> / dq 2 \q2){n\ , 
JS 2 



Xi = / dp 1 \pi)(pi\ <g> 1 2 , 

X2 = ll ® / dp 2 \P2)(P2 \ ■ 



(IV.5) 



The x's are orthogonal projectors^ (xt = X> X? 2 = x) acting on TL = L 2 (M.,dq\) ® 
L 2 (R, (iq^)- The product of two of them involving different indices commutes, so that 
V is a (bounded) self-adjoint operator. 

The inequalities (III. 2) to be tested in the quantum context then become, for pure 
states p = |^)(*|, 

0<(*|P|¥)<1 V|*)eWwith <*|*) = 1, (IV.6) 



3 In eqs.(EV.5), Si, S2, S[, S' 2 are the supports of xi> X2, Xii X2 respectively. Also, 
J s dqi \qi){qi\ is, in standard Dirac notation, the orthogonal projection ip(qi) — ► Xi(9i)V'( < /i)! 
whereas J" g , dpi |pi)(pi is the orthogonal projection ip(pi) — * Xiip^'ipiPl)) $(Pi) being the Fourier 
transform of tp(qi), and so on. 
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or, equivalently, 

V > and 1 — "P > in the operator sense. (IV. 7) 

Because Xj f aus to commute with x'j (j = 1)2), "P is not an orthogonal projector 
(see below), in contrast to the classical equality V 2 = V. Exploiting this fact leads 
to the 

Proposition 1 The operators V and (l — V) cannot be both positive. 

As a consequence, there is at least one \^>) ^ such that the inequalities > 

and (^1(1 — V)^} > cannot be simultaneously true. This just means that one of 
the two inequalities (IV.6) is violated for that \^), which settles the question. 

Proof of proposition 1: 

Assume that V and (1 — V) are both positive. This would imply 

V(1-V)>0, (IV.8) 

(remember that the product of two positive commuting operators is positive). 
Now, a straightforward calculation of V 2 from eq.( pV.4 ) yields 

V 2 = V- [xi,xi] [X2,X 2 ] , (IV.9) 



and eq.( IV.8 ) would mean that [xi>x'i] [X2>X 2 ] i s a positive operator. That this is 
wrong is not surprising. Let us show it. Take a factorized \^>), namely |\&) = |3>i) <g> 
| $2), so that 

(y\V(l-V)\y) = -($ l \i[x 1 ,X , i] |$1)($2K[X2,X 2 ] 1*2). 

It is enough to show that, for a given choice of the characteristic functions x an d x! ■> 
the real number 

R[$\ = <*|i[x,2l 1$) (IV.10) 
can assume both signs when |$) is varied. 
Let us define 

i$+) = xi$), 1$-) = (i-x)m- 

Using the identity 

[x, x] = xx'i 1 - x) - (1 - x)xx 

gives R\&\ the form 

R[$} = i($+|x'|$~) - i(^~\x'\^ + ) • 
Obviously, for |$) = |<I> + ) - l^ - ), one has R[$] = —R[$]. 

This concludes the proof. 
Remarks: 

1) When the wave function factorizes, i.e. ^(<7i, (72) = $i(<7i)$2(<72)j a correspond- 
ing probability distribution p always exists, namely 

p(<7l,<B,Pi,P2) = |$i((?i)| 2 |$ 2 fe)| 2 |$i(R)| 2 |^2(p2)| 2 , (IV.ll) 
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where the $j's are the Fourier transforms 

1 r+oo 



(IV.12) 



/2ir j-oo 

Of course, this implies that eqs.(IV.6) are automatically satisfied for such factorized 



|^)'s (which can also be checked from eq.( lV.4 )). 

2) The fact (used in the proof) that (^^(l — < for some factorized |^)'s is 

not inconsistent with the inequalities < ( 1 I'|'P| , I / ) < 1 which are satisfied for those 
|*)'s. 

IV-B. Construction 



We want to find wave functions |^) violating the inequalities ( IV.6 ). According to 
the first of the above remarks, one has to depart from the class of factorized |\&)'s. 
The simplest way to do it is to take just a sum of two such products. 

Choose first 

S\ = S2 = S and S[ = S' 2 = S' , 



so that 



^ = X®l2 + li®X + X ® x - x®x- x®x -x ® X ■ 



(IV.13) 



Take next 



*) = ^i^(|0) + A^)) (AeC), 

|0) = |0i> ® |0 2 ) , \<p) = \<pi) ® \<f>2) , 

(01 1 01 ) = (02|02) = (<Pl\<Pl) = (^2^2) = 1, 



with 



bi) =0, 



so that is properly normalized. 
For the moment, choose also 



with 
Then 

with 



1 + A 



1 = 02 = / and ipi = ip 2 = g ■ 
(f\f) = (g\g) = l, (f\g) = 0. 

t>\V\<t>) + (A(0|P|^ + c.c.) + \\\ 2 (<p\f>\ip) 



(IV.14) 

(IV.15) 
(IV. 16) 

(IV.17) 



(01%) = 2(/| X |/) + (/| X '|/) 2 -(/| X |/) 2 -2(/|x|/)(/|x'|/) 
{<p\P\<p) = 2(g\ X \g) + (g\ X >\g)Z-(g\ X \g) 2 -2(g\ X \g)(g\ X >\g), 

(<t>\V\v) = (f\x'\g) 2 -(f\x\g) 2 -2(f\x\g)(f\x'\g)- 



(IV. 18) 
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We already know that < {(f)\V\(t)) < 1 and < (ip\V\(p) < 1. Clearly, in view of 
(IV. 17), our goal will be reached (namely (^"Pl^) < or (\E'|'P| V I / ) > 1) if one can 
find / and g such that 



(IV. 19) 



We claim that this can be achieved with S = S' = (0, oo), f(q) = (q\f) an even, 
normalized function in L 2 (— oo,oo), and 



g(q) = (q\g) = sgn(g)/(g) . 

With this choice, eqs.( IV.l^ ) are automatically satisfied and 

(f\x\f) = (g\x\g) = (f\x\g) = l- 



(IV.20) 



(IV.21) 



Also, since the Fourier transforms f(p) and g(p) are respectively even and odd func- 
tions 

(f\x'\f) = (9\x\g) = l- (IV.22) 

As for the non trivial interference term (/|x'|s)> it is given by (see appendix A) 

P 



(f\x'\g) 



7T JO 



OO POO f 1 



q + q' q — q 



(IV.23) 



At this stage, it is advantageous to take / as a reai function, so that by symmetry 

, mm 



(f\x'\g) 



7T Jq 



OO ("OO 



dq \ dq 



Let us set 



Then 



with 



and 
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h(q) = V2f(q), 
K(q,q') = 

vr q + q' 

(f\x'\g)=-^, (7€R) 

roo 

dq / dq'h(q)K(q,q')h(q') 



lL 2 (0,oo) — 1 • 

The insertion of eqs. ([lVT2"l"l) , (|rV2l ) and ( |TV26l) in (fVT8D gives 

(<f>\P\<t>) = (<p\V\<p) = l, 
(0|P|^) = -i(l + 7 2 ) + |7, 



(IV.24) 
(IV.25) 

(IV.26) 
(IV.27) 
(IV.28) 



so that eq. (|IV.19D reads 



( 7 2 + l) 2 + 4 7 2 >4, 
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which is satisfied provided that 



|7| > y2V3 - 3 0.6813 
Moreover, with X = pe ie , eq.( IV.17| ) becomes 

1 P 



2 2(1 + p 2 



(1 +Y) cos # + 27 sine 



(IV.29) 



(IV.30) 



We already know that I7I cannot exceed 1, because |(/|x'|#}| 2 < (f\x'\f) (g\x'\g) = \- 
Are there however some h's (subjected to ( IV.28] )) such that 7 (given by eq.( [IV.27| )) 
fulfils ( IV. 29 )? If this occurs we have reached our goal and it only remains to maximize 
I7I in order to obtain the extremal values of ( 1 I / |'P| 1 I / ) (within the present scheme) 
through eq.( IV.30| ). In other words, one has to solve the problem 



7o 



sup 

L 2 (0,o< 

h=h* 



\(h\K\h)\ 



'l 2 (0,oo) 1 



In appendix B, it is shown that the (bounded) integral operator K with kernel 
( IV.25| ) on L 2 (0,oo) is positive and has the purely continuous spectrum [0,1]. This 
immediately entails 70 = 1, and we get 



7=1 



t (cos + sin ( 

1 + p z 



inf mv\$) 

A 

sup {y\V\% 

A 



7=1 



7=1 



7=P=l,9=f 



7=p=l,6 



-3tt 
4 



2 

1+V2 
2 



-0.2071 
1.2071 



(IV.31) 



Actually, as discussed in appendix B, due to the continuous spectrum of K, these 
extremal values cannot be strictly reached, but only approached arbitrarily close via 
a family of normalized functions h, e.g. 



M?) 



9{L 



1 



v/ln(L + l) V^TT' 



L 



DC 



or smoothed forms of this. Of course, other functions h will also do the job (although 
less perfectly), that is meet the crucial requirement ( IV.29| ). Taking for example 
M?) = ^XT (which is normalized in L 2 (0,oo)), one gets 



T = 0.7854 
4 



Finally, collecting the equations ([TV.14P , (|IV.15| ), < \W.2(\ ) and (jTV.240 , togeth er with 
A = f ; — X' one obtains the wave functions leading to the maximal violations ( IV.31 ) 



2^2 



1 ± e*4 sgn(gi) sgn(g 2 ) /i(|9i|)/i(|9 2 |) 



(IV.32) 



where h(q) stands for some regularized form of with f£° dqh(q) 2 = 1. 
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IV-C. Introducing Einstein locality and relative motion 



Let us now interpret q\ and q 2 as the coordinates of two particles (rather than the x 
and y coordinates of the same particle).. Then the wave functions (|IV.32| ) describe 
states of two particles not spatially separated and with zero relative momentum. 
These two restrictions can be easily disposed of. 

First, it can be checked that nothing is essentially changed in the previous deriva- 
tion if one keeps 

51 = £l = (0,oo); = Pi(gi) = sgn(gi)/(gi), (IV.33) 
but replaces 

5 2 = S 2 = (0, oo) ; <f> 2 (q 2 ) = /fe) , ^2(92) = sgn(q 2 )f(q 2 ) , (IV.34) 

by 

S 2 = (a,oo), S' 2 = (0,oo); (f) 2 (q 2 ) = f(q 2 - a ) , (f 2 {q 2 ) = sgn(q 2 - a)f(q 2 - a) . 
Then eq.( iV.32] ) becomes 



2V2 



l±e l i sgn(gi)sgn(g 2 - a) MM)MI<72 - a\) 



with a arbitrary. 

This allows us to let Einstein locality enter the game. 

Similarly, nothing is essentially changed if one keeps eqs.( pV.3"3 ) but replaces 
eqs. (HV.34j) by 



S 2 = (0, 00) , S' 2 = (P, oo) ; MQ2) = e iP<?2 /fe) , ^2(92) = e iP " 2 sgn(g 2 )/(g 2 ) . 
Then eq.( IV.3^ ) becomes 



*± (91,92) 



1 ± e*i sgnfe) sgn(g 2 )l e iP< *h(\ qi \)h{\q 2 \) , 



with P arbitrary. 

This allows us to put the two particles in relative motion. 

IV-D. Analogy with Bell spin | correlation inequalities 

Let us denote by |+) a normalized function / close to the (symmetrized) eigenfunc- 
tion of the operator K with "eigenvalue" Ao = 1 (i.e. 7 = 1 in eq.( [IV.26 )), and 
by |— ) the orthogonal function g (as given by eq. (|IV.20|) ). Consider the subspace 
V = span(|+), |— )) of the full 1-particle Hilbert space, together with the orthogo- 
nal projector II onto V. Call T (resp. V) the restriction of x (resp. x') to the 
2-dimensional space V 

r = nxn, r' = nx'n. 
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Then eqs.( |[V.2lD , flTV.22|) and (|W.26| ) tell us that T and T' are represented in the 
orthonormal basis {|+), |— )} by the matrices 



r 



r 



1 i 

2 2 

1 1 

2 2 

37 



3 (! + **)> 



\ (1 - 7CT y ) . 



In the idealized limit 7 — ► 1 (and only in this limit), one observes that T and V are 
themselves orthogonal projections V — > V 



r 
r' 



rt, 
r't 



r 2 
p/2 



r, 
r. 



This implies that both operators x an d x' leave the subspace V invariant 

[n,x] = [n,x'] = o. 

Indeed, a straightforward calculation shows that 

[(i-n)xn]t[(i-n)£n] = o, 

which entails (1 — II)xII = and xll = Tlx- The same for x' ■ 

Hence, in the 2-particle Hilbert space, the operator ( |IV.13 ) also leaves invariant 
V <S> V, and V := "PII assumes the simple form 

(Vl) ff (2) _ a W a m + CT (l) a (2) + a (l) a (2) 

\ y y x x x y y x 
whereas the maximally violating wave functions ( [IV.32 ) read 

, : - — M+V l '|-i- •)'•-' +,.'Tl-\( 1 )|-\( 2 )' 



- 1 1 

V = - + - 

2 4 



(IV.35) 



-^(|+)«|+)( 2 )±e^| 



From eq.(!V.35) one can check that 

V(l-V) 



(IV.36) 



(IV.37) 



which is just the projected form of V(l — V) = [xi, Xi] [X2> X2]' an< ^ the expectation 
value of the operator ( IV.37|) is 

(*±F(i-P)|*±) = -i 

for the wave functions ( [IV.36D . 

The result ( IV.31| ) is also directly recovered from eqs.( |IV.35 ) and ( IV.36| ) 

Then one sees that, in the idealized limit 7 — * 1, the original phase space setting 
up of the problem is formally equivalent to the standard EPR setting up for a two 
spin I system, together with its classical Bell inequalities. 
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V. General solution of the three marginal problem 



We have proved here the impossibility of reproducing quantum probabilities of four 
CCS as marginals. Roy and Singh have given examples to show that reproduc- 
ing three CCS is possible. In this section, we construct the most general nonnega- 
tive phase space density which reproduces three different (noncommuting) CCS as 
marginals. Our results encapsulate the extent to which noncommuting CCS can be 
simultaneously realized in quantum mechanics. 

Among the four marginals R, S, T, U obeying the compatibility conditions (|II.7| ) 
which are at our disposal, the particular choice of three of them is completely irrele- 
vant. For defmiteness, we choose R, T and U, which we rename 00(91,92), 0-1(^1,92) 
and a 2 (pi,P2)- 

We assume that these marginals are probability densities in the full mathematical 
sense, that is they are true (integrable and non negative) functions. This means that 
we restrict our marginal probability distributions to absolutely continuous measures 
(with respect to Lebesgue measure) in R 2 . Notice that such a restriction is automatic 
in the quantum case, due to eqs.(II 



Likewise, we look for the general solution of the three marginal problem in the 
class of absolutely continuous measures in the phase space R 4 . This means that we 
want to describe all the solutions p of the equations 



co(<7i, 92) = J dp 1 dp 2 p(~ct, ~f) , 

92) = J dqidp 2 p(lt, y) , (V.l) 
0"2(Z>1,P2) = J dqidq 2 p{ if, , 



which belong to L 1 ^ 4 , d 2 qd 2 p). 

Notice that this is a restricted problem even in the quantum case, since nothing 
prevents a probability measure containing a singular part to project on marginals 
which are L 1 -functions. To some extent, the above restrictions can be removed, al- 
lowing us to include e.g. probability measures partly concentrated on submanifolds 
of the phase space. However, dealing with such extensions at some degree of gen- 
erality requires painful manipulations, and we shall ignore them here0. As for the 
full inclusion of singular measures, it appears as both delicate and of little practical 
interest. 

Let us introduce the one variable marginals 



0-01(92) = / dqi 0-0(91,92) , 

K (V.2) 
fl2 (Pi) = / dq 2 (7 1 (pi,q 2 ) . 



'Special cases are treated in M and 
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Owing to the compatibility conditions (|IL7 ), these definitions are equivalent to 



ooife) = J dp! <ri(pi, q 2 ) , 

Vnipi) = J dp 2 (T 2 {pi,P2) ■ 



(V.3) 



As the support properties of the functions aj (which are allowed to vanish on 
some parts of R 2 ) are not innocent in the forthcoming construction, we need to pay 
attention to them. Let Ej C R 2 (j = 0, 1, 2) be the essential support of Oj. The above 
compatibility conditions, together with the positivity conditions aj > 0, clearly yield 
two constraints on the supports Ej, namely 



(V.4) 



{q 2 G R | 3 q 1 G M such that (qi,q 2 ) G E } = 

{q 2 G M | 3 pi G M such that (pi,q 2 ) G Ei} , 

{pi G R | 3 q 2 G R such that (pi,q 2 ) G Ei} = 

{pi G R | 3 p 2 G R such that (pi,_p 2 ) G E 2 } . 

To the Ej 's we associate the subsets Ej 's of the phase space defined by 

Eo = {?, f | (gi,g 2 ) G E ,(pi,p 2 ) G R 2 }, 

^i = {"?, I (Pi , 92) G E x , ( gi , pa) G R 2 } , (V.5) 

E 2 = f I (pi,P2) G E 2 ,(q h q 2 ) G R 2 } . 

Finally, we denote by E the intersection of the Ej 's 

£ = £ n E 1 n E 2 . (v.6) 



Clearly, due to positivity again, any solution p of eqs. (|V.l|) must have its essential 
support contained in E. 

The three marginal problem in the precise form stated above is then completely 
solved by 

Theorem 11) The Lebesgue measure of E is not zero and the function po defined 
(a.e.) by 

f O"o(gl)g2) - ) n y cri0l,g2) - ) n v (T 2 (pi,p 2 ) if (~ct,l?) G E , 
[ otherwise, 

is a non negative solution of the problem ( \V-1\ ) in L 1 (M. 4 ,d 2 qd 2 p). 
2) The general solution p of in L 1 (R 4 ,d 2 qd 2 p) is given by 

Pit, f) = Po(t, P > ) + AA(t, f), (V.8) 

where 



A G 



1 1 

m + ' m_ 



(V.9) 
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and 

A(~<f , ft) = F(lt, ft) - po(!t, P~) 
+ 



1 



o"l(Pl,92) 
1 

0-01(92) 



0-0(91,92 
dq'tdp'z F(q[ , q 2 , pi, p' 2 ) + 

dq\dp\dp 2 F{q[ , q 2 ,p[ , p' 2 ) - 



dp' 1 dp' 2 F(q 1 ,q 2 ,p' 1 ,p 2 
1 



02(^1,^2) 
1 

012 (pi 



d4xdJ 2 F{qU*,PX,pi) 
dqxdq 2 dp 2 F(q' 1 ,q 2 ,p 1 ,p 2 ] 



(v.io) 

F being an arbitrary L 1 (M 4 ,d 2 qd 2 p)- function with essential support contained in E. 
The (F- dependent) constants m± in ^V.dJ are defined as 



m+ = ess sup — 

(T,-?)£E Po{Y , P ) 



. , A( If, ft] 
m- = — ess mj — _^ 

("g > ,"P > )6S M 9 > P 



(V.ll) 



and are k>i/i positive if A 7^ (m + = 00 or/and m- = 00 are not excluded). 
Proof: 



1) To begin with, po given by ( V.7 ) is well defined and non negative. Indeed, due to 
( |V. 2D (or (|y.3|) ) and the positivity of the cry's, o-Q\{q 2 ) and ai 2 (j>i) are a.e. non zero 
for (if , ft) £ Eq and £Ji (or E\ and ^2), so that the denominators in eq.( V.7 ) do not 
vanish on E (except maybe on sets of Lebesgue measure 0) . 

Next, in order to check that po obeys the first equation ( V.l ), we consider the 
integral 

dpi / dp 2 p (lt,ft) 



(V.12) 

with this specific order of the p integrations. According to the relations (|V.4| ) and 
the definition of E, one observes first that the projection of E on the (pi,p 2 ) plane 
is the set Y, 2 , so that the integration over p 2 removes the factor a 2 ja\ 2 in po! an d 
second, that the projections of Si and T, 2 on p\ coincide, so that the integration over 
p\ removes the factor cji/cjoi in po, and one is left with the expected result a"o((?i, q 2 ). 
We can now write 

dp 1 dp 2 p Q (~q*,ft) = 0-0(91,92) (V.13) 



where, thanks to Fubini theorem, the integration order is completely irrelevant. The 
other two equations (V.l) are derived in a similar way. 

This calculation shows at once that the Lebesgue measure of E is not zero and 
that po G L l (R 4 ,d 2 qd 2 p). 

2) That any non negative solution p of eqs.( |V.l| ) admits the representation ( |V.S| )- 
flV.lCD is easy to establish. Indeed, since the essential support of p is necessarily 
contained in E, we are allowed to take F = p in eq.( |V.10| ), which gives (using (V.l)) 



A(lf , ft) = p{~t , ~f) ~ Po( It, ~tf 



(V.14) 



Then, from (|VT1 ) 



m_ 



ess inf ( — 



1 <1. 
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As l/m_ > 1 in (|vTg|) , we can choose A = 1, which makes eq. (|V.14| ) equivalent to the 
representation ( |V.8| ). 

It remains to show that any function p defined by (|V.8| ) to ([V.ll] ) (and thus with 
essential support E) is a non negative solution of eqs. (fVTT]) in L 1 (R 4 , d 2 qd 2 p). In order 



to prove that p satisfies the first equation ( V.l ), we rearrange pairwise the right-hand 
side of ( |V.10 ) as follows 



A 



F 



dp[dp' 2 F 



d4\dp\dp' 2 F 



_<T\ J 001 

— / dq[dq 2 F- — / dq\dq 2 dp 2 F 

.0"2 J <7i2 J 

Then, integrating the right-hand side over p\ and P2, one finds, by an extensive use of 



(V.15) 



eqs.( V.2 ) to ( V.4 ) as in part 1), that the two terms coming from each square bracket 
cancel each other, leading to 

dpidp 2 A(~q*,J?) = 0. 



This, with ( |V.8 ) and ( V,13| ), implies that p satisfies the first equation (V.l). That it 
satisfies the other two equations ( V.l ) is proved in a similar way. 

This calculation also shows that p € L 1 (R 4 , d 2 q d 2 p). 

Finally 

/ d 2 qd 2 pA(lf,jf) = 0, 
JE 

which implies that m± in eqs. QV.il ) are both strictly positive if A does not vanish 
a.e. on M 4 . The positivity of p is then a trivial consequence of eqs.( |V.8| ), ( |V.S| ) and 

The proof is complete. 
Remark: 

Theorem 1, as it is stated above, deals with L 1 functions, and thus excludes the 
occurrence of Dirac measures. We insist on the fact that this is unnecessarily restric- 
tive. Indeed Dirac measures can be easily accommodated and the theorem suitably 
rephrased, to the price however of cumbersome mathematical intricacies which we do 
not want to enter into. 

An immediate corollary of Proposition |l| and Theorem [l] is 

Theorem 2 (Three marginal theorem) Let R, S, T and U be probability distri- 
butions for (q\,q2), (qi,P2), (pi,Q2) and (pi,P2) obeying the consistency conditions 



(H.7). Given n arbitrary distributions among {R, S,T, U}, a necessary and sufficient 
condition for them to be marginals of a probability density in the 4- dimensional phase 
space is n < 3. 



VI. Conclusions 



We have solved the four marginal problem in four dimensional phase space thus 
proving a long standing conjecture P| and vastly improving the first results of Martin 
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and Roy which dealt with infinite number of marginals. To achieve this, we first 
derived "phase space Bell inequalities" which have their own interest. Actually they 
allow, at least in principle, direct "experimental" tests of the orthodox-versus-hidden 
variable interpretations of quantum mechanics within the position-momentum sector, 
analoguous to those performed within the spin sector. 

The technique of phase space Bell inequalities established here has applications 



to quantum information processing. Generalizing the example ( IV. 11 ), one can show 



that for any separable density operator p one can construct a phase space density 



obeying the four marginal conditions. Hence, the Bell inequalities ( [II. 2 ), with R, S, 
T and U given by ( |II.8| ) must hold for every separable quantum state, irrespective of 
any physical interpretation of the associated phase space density. Their violation by 
a quantum state is a signature and even a quantitative measure of entanglement of 
this state. 

We have also constructed the most general positive definite phase space density 
which has the maximum number of marginals (three) coinciding with corresponding 
quantum probabilities of three different (noncommuting) CCS. These results should 
be useful in the construction of maximally realistic quantum theories. 



VII. Acknowledgements 

We thank Andre Martin for collaboration in the initial stages of this work. One of us 
(SMR) thanks A. Fine and A. Garg for some remarks on the three marginal problem 
many years ago. 



17 



Appendix A. Proof of equation fliV 

Since S' = (0, oo), one has 

x'{p)g(p) = 8{p)g{p) ■ 

Assuming first that g belongs to S (the Schwartz space of infinitely differentiable 
functions on M with fast decrease at infinity), one can write 

(x'g){q) = / dq 9{q - q')g(q') , 



where is the Fourier transform of 9 in the distribution-theoretic sense 
%) = oz/ *e«0(p) = -- + -%). 

27T 7_oo 27T q 2 

Then, if / also belongs to S 

(/|X b) = f / d?/* (?) / dq' 7 g(q') + - (f\g) . 

2tt J~oo J-oo q — q 2 

In particular, for even / and odd g, (f\g) vanishes, and 

i r°° r°° ( 1 P \ 
{f\x\g)=— dqf{q) dq' ( —— ) g(q') , 

7T Jo Jo \q + q q-qJ 



which gives eq.( IV.23| ) if g coincides with / on (0, oo). The continuation from S to 
L 2 (— oo,oo) is performed as usual by continuity, using the fact that S is a dense 
subspace in L 2 ( — oo,oo). 



Appendix B. Study of the operator K 



From the very definition of K through the integral kernel ( IV.25| ), one has 

n Jo q + q 

Let us put 

h(u) =e%h(e u ). 

Since /f^ du \h(u)\ 2 = J °° dq \h(q)\ 2 , the correspondence h i— > h defines a unitary 
mapping L 2 (0, oo) — ► L 2 (— oo,oo) and 

/oo 
du if(it - v)h(v) , (B.l) 
-oo 

where 

K(u) = . 

y ' 2vrcoshf 

Then, another unitary mapping L 2 (— oo,oo) — > L 2 (— oo,oo), namely the Fourier 
transform 

/i(fe) = -= / due iku h(u), 
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reduces the convolution product in (|B.l| ) to an ordinary product 

Kh(k) =K(k)%(k), 

where 

_ /■oo 1 

K(k) = / due lku K(u) = — — - . (B.2) 
J-oo cosh7r/c 

Therefore, the operator K on L 2 (0, oo) is unitarily equivalent to the multiplicative 
operator ( |£?.2| ) on L 2 (—oo, oo). The latter is evidently a positive operator with purely 
continuous spectrum [0,1]. Its generalized (non normalizable) "eigenfunctions" are 

h s (k) = 5{k - s) (s£l), 
with "eigenvalues" X s = — ^ ^ , and their preimage in L 2 (0,oo) are 

Of particular interest for us is the extremal one, with "eigenvalue" Ao = 1 

ho{q) 



Of course, the corresponding maximal value 70 = 1 of 7 = (h\K\h) cannot be attained, 
but only approached arbitrarily close through a family of normalizable functions mim- 
icking -^=. For instance, introducing two cutoffs, e at small q and L at large q, and 
setting 

he,L(q) = —j= X{e,L)(q) ~p = 1) , 

Jin f V « 

one gets 

M k Kl ) = 1 - 2^ = 1 - o ( ' ) , 



7rln-^ Jv/i7^ x \ln 



e v ' \ e / 

so that lim (/i £ i\K\h £ l) = 1 • 

L^oo 

Notice that one can keep e fixed (e.g. e = 1) and let L alone go to cxd without 
changing anything (this is in fact a consequence of the scale invariance of the operator 
K), or even choose a family of less singular functions h, like 

h L {q) = -=±==9(L 



y/HL + 1) vWT' 
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